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ABSTRACT 

We compare broad emission line profiles and estimate line ratios for all major emission 
lines between Lya and H/3 in a sample of six quasars. The sources were chosen with 
two criteria in mind: the existence of high quality optical and UV spectra as well as 
the possibility to sample the spectroscopic diversity in the 4D Eigenvector 1 context 
(4DE1). In the latter sense each source occupies a region (bin) in the FWHM(H/3) 
vs. Fe iiopt strength plane that is significantly different from the others. High S/N 
H/3 emission line profiles are used as templates for modeling the other lines (Lya, 
CivA1549, HeiiA1640, AhiiA1860, Siiii]A1892, and MgiiA2800). We can adequately 
model all broad lines assuming the existence of three components distinguished by 
blueshifted, unshifted and redshifted centroids (indicated as blue, broad and very 
broad component respectively). The broad component (electron density Uc^ 10^^ 
cm~^; ionization parameter U ^ 10~^; column density Nc> 10^^ cm~^) is present 
in almost all typc-1 quasars and therefore corresponds most closely to the classical 
broad line emitting region (the reverberating component). The bulk of MgiiA2800 and 
Fell emission also arises in this region. The blue component emission (log ric^ 10; 
log (7 ~ — 1; log iVc< 23) arises in less optically thick gas; it is often thought to arise 
in an accretion disk wind. The least understood component involves the very broad 
component (high ionization and large column density) which is found in no more 
than half (but almost all radio-loud) type-1 quasars and luminous Seyfert nuclei. It 
is perhaps the most distinguishing characteristic of quasars with FWHM H/3>4000 
km s~^ that belong to the so-called Population B of our 4DE1 space. Population A 
quasars (FWHM H/3<4000 km s~^) are dominated by broad component emission in 
H/3 and blue component emission in CivA1549 and other high ionization lines. 4DE1 
appears to be the most useful current context for revealing and unifying spectral 
diversity in type-1 quasars. 

Key words: galaxies: quasars - galaxies: quasars: general - galaxies: Seyfert - line: 
profiles - line:formation - galaxies: quasars: individual: I Zw 1, Mark 478, Mark 335, 
Fairall 9, 3C 249.1, 3C 110 



1 INTRODUCTION 

1.1 Introduction 

Quasar spectra are not self-similar even if we restrict our- 
selves to active galactic nuclei that show broad emission lines 
(type-1 AGN). Examination of the same broad emission line, 
in many different sources reveals a striking diversity. This di- 
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versity begins with the profile width but extends to every 
other measurable property (e.g. shape, line shift, equivalent 
width). The Hi Balmer line H/3 is a good place to start a 
line profile characterization because so much data exist. It 
can be observed from the ground out to z;~0.9 and through 
numerous IR windows to z >3. Thus one can compare the 
line over a wide range of redshift and source luminosity. It 
is usually a strong line and suffers only moderate contam- 
ination falling in a gap between strong optical Fen blends. 
At this time good S/N and moderate resolution spectra ex- 
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ist for more than 600 low ^ sources (|Marziani et al.lll99q . 
l2003bl : IZamfir et al. 200 8, 201^. includinR t he brig htest «400 
SDSS quasars). Table 2 in ISulentic et all (|2000al ') tabulates 
earlier references with high S/N spectra. Comparable qual- 
ity IR spectra of the H/3 region n ow exist for more than 50 
sources in the range z = 0.9 - 3.1 (|Sulentic et al.ll2004[200^ : 
iMarziani et al. 2009 ). 

The Civ line at H549 is the other current broad line- 
of-choice providing a high ionization (HIL) counterpoint to 
low ionization (LIL) H/3. High quality ClvA1549 spectra for 
low z sources are more difficult to obtain with the HST 
archive provi ding UV spectra of the ClvA1549 region for 
~130 sources (JBachev et al.ll2004l '). Above z ~1.5 CivA1549 
is accessible at optical wavelengths out to z ~4.8. The 
largest compilation of high quality and high z ClvA1549 
spectra came from the Palomar surveys almost 20 years ag o 



l| Young et al.ll 19821 : ISargent et aLlll989l : iBarthel et a"Llll990 
A uniform study of the Palomar and other g ood quality HIL 
spectra can be found in lTvtler fc FanI (|l992l '). Unfortunately 
most sources with good H/3 spectra lack a corresponding 
ClvA1549 spectrum and vice versa. Many of the other broad 
lines have been studied but usually without the benefit of 
spectra for other broad lines in the same source. A funda- 
mental desideratum involves strong narrow lines that pro- 
vide a reliable determination of the quasar rest frame (e.g. 
[Oiii]AA4959,5007, and especially [Oii]A3727). Full exploita- 
tion of the Palomar spectra is impeded by the lack of rest 
frame determinations for the sources and the lack of match- 
ing H/? spectra. 

Broad line profiles provide some of the most direct 
clues about the nature, structure and kinematics of the cen- 
tral regions in quasars. Inter- line profile comparisons for 
individual sources (Seyfert 1 nuclei as well as radio-quiet 
(RQ) and radio-loud (RL) quasars) have been carried out 
for a long time alt hough there has been litt l e agreement 
on the systematics (lGaskelllll982l:IWilkeslll98llEspev etla] 
1989':'C orbinlll99Gl :'Carswcll et al."l99l':'T vtler fc Fanlll992r 
Marzia ni et al.l ll996: Richards ct al. 200^. Most of these 
studies have sampled: (a) small numbers of sources, (b) 
only a few of the broad l ines in each source, and/or (c) less 
than ideal spectra (see IShang. Wills. Wills, fc BrothertorJ 
((20071) for a notable exception). If we consider the best 
studied HIL and LIL features (ClvA1549 and H/3) we 
find little evidence for profile (shape/shift/EW) similarity 
for sources with FWHMH/3< 4000 km s~^ while above 
that value the profiles appear to be more similar, af- 
ter proper corre ction for narr o w line ClvA1549 emission 
made (iMarziani et al.l 1 19961: ISulentic fc Marzianil 1 19991 : 



iBachev et al.ll2004l : ISulentic et al.l 120071 '). Ideallv we would 
like to have spectral coverage for all of the principal broad 
lines in each of a sample of sources that reflect the impres- 
sive source-to-source diversity now known. It would also be 
ideal if we could achieve a consensus on how many distinct 
components each line contains. This paper will address both 
of these desiderata. 

1.2 Eigenvectors of Quasars 

We are never satisfied with the S/N and resolution of the 
available quasar spectra even though current best data ex- 
ceed our ability to explain them theoretically. Minimum ac- 
ceptable values for good spectra are arguably 20-30 in the 



continuum and A/AA ~ 1000 respectively. We are tempted 
to indiscriminately average line profiles in order to improve 
even more the S/N but a composite spectrum generated 
from a large sample obscures most of the important astro- 
physical information. Such composite spectra have heuristic 
value since they provide a wide spectra l coverage and give 
a glo bal view of the continuum shape (jVanden Berk et al.l 
[200l|) . Different lines represent averages of different sources 
making comparison impossible unless all quasar spectra are 
fundamentally similar. In the context of broad line physics 
indiscriminate averaging obscures rather than reveals profile 
differences. It is the "why" of differences that yields physical 
insight. Not only do quasar spectra show significant diver- 
sity but measured line parameters do not scatter symmet- 
rically around an average. Eigenvector analysis showed that 
the diversity of optical and UV spectra i s not entirely ran- 
dom, but that there are syst ematic trends (jBoroson fc GreenI 
1 19921 : ISulentic et al.|[2000bh . We expanded the original opti- 
cally based Eigenvector 1 parameters to i nclu de ClvA1549 
and soft X-ray photon index (see also iWang et al., il996l ) 
in an effort to enhance our ability to identify profile d if- 
ferences (4D Eigenvector 1 = 4DE1 ISulentic et al.|[2000al lbl. 
[2007]) . As examples of the lack of symmetric scatter we note 
that sources occupy well-defined non-random sequences in 
e.g. FWHM(H^) vs. normalized FenA4570 strength (see Fig- 
ure [T]) or vs. norma lized ClvA1549 line shift (Figure 2d in 
ISulentic et"alll2007l ). 

Profile diversity persists when composite spectra are 
generated in some fundamental context like 4DE1. For ex- 
ample, sources in a restricted region of the 2D parameter 
space involving FWHM H/3 and Fe llopt strength measures 
show considerable profile diversity (Figure [T] which is the 
optical plane of 4DE1). The composite spectrum for sources 
in a restricted region (bin) presumably reveals the under- 
lying "stable" emission component (s). The source-to-source 
diversity within a bin likely refiects minor changes in pro- 
file shape connected with profile variability. The bin-to-bin 
diversity in a contextual space like 4DE1 have the great- 
est potential to yield clues about broad line region (BLR) 
structure and kinematics (convolved with source orienta- 
tion). Contexts for such binned profile composites can be 
empirical (e.g. FWHM, Fen strength) or physical (e.g. black 
hole mass Mbh, Eddington ratio I/boi/iEdd). A comparison 
between e.g. radio-quiet and radio-loud quasars will yield 
any result that one wants unless the empirical context for 
profile binning is carefully chosen and is independent of ra- 
dio properties. 



1.3 How Many Broad Line Components? 

We would like to compare all principal broad lines in a small 
sample of sources that reflect the diversity found in the 4DE1 
binning. At present however adequate multiwavelength data 
exist for very few sources so studies encompassing all the 
principal broad lines are rare (e.g. Shang et al. 2007). The 
most basic question for facilitating such a study involves how 
many components one can identify in each broad line. This 
is particularly important for studies that use the widths of 
broad lines as virial estimators. Empirical studies have made 
it clear that HIL (e.g. ClvA1549) and LIL (e.g. Balmer se- 
ries) show a wide range of strength, shape and shift prop- 
erties. Our own work in the 4DE1 context has attempted 
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to systematize these properties iSule ntic et al.ll2000al . 120071 : 
iMarziani et~ai1 1200 ll . |2003b c. 2008) with focus so far on 
comparisons of ClvA1549 and H/3 assumed to be typical 
HIL and LIL features. 

Several authors have considered the degree of strati- 
fication that might exist in the BLR including the idea 
that the innermost part of the BLR could be emitted in 
a quasi-distinct very broad line region (VBLR). The VBLR 
concept may have originated with the discovery that the 
core and wings of HIL HellA4686 and LIL II/3 in highly 
reverberated NGC5548 show different re sponses to contin- 
uum change. (|Peterson fc FerlandlllQSeD . The VBLR con- 
cept and/or terminology arose in several additional guises 
and this has led to some confusion. Is the VBC present 
in all quasars? Are t here equal nu mbers of blue- and red- 
shifted VBC (see e.g. ICorbinlll995l )? Is there no significant 
narrow emission line component in ClvA1549? If one rejects 
a significant ClvA1549 narrow component the "peakiness" 
of the ClvA1549 profile in many sources suggests both nar- 
row intermediate line region (ILR) and very broad VBLR 
as th e simplest ClvA1549 profile model (jBrotherton et al.l 
1 19941 '). W e think the answer to a ll of these questions is "no" 
(see e.g. IMarziani et al.l 1 19961 : ISulentic fc Marzianil 1 19991 : 
ISulentic et al.ll2002l . |2007^ . Our studies show evidence for a 
VBC only i n sources with FWHM (H/ 3) > 4000 km s~^ (Pop- 
ulation B: ISulentic et al.l l2000al rbl: IMarziani et all l2003d : 



IZamfir et al.ll2010l ). The II/3 profile in these sources requires 
both unshifted gaussian (BC) and redshifted VBC compo- 
nents. Narrower (FWHM II/3< 4000 km s"^) Pop. A sources 
can often be fitted with a simple Lorentz-like function al- 
though a blue asymmetry is seen in the H/3 profiles of some 
pop. A sources. Whether FWHM H/3« 4000 km s"^ signifies 
a fundamental change in BLR physics (e.g. critical Edding- 
ton ratio?) or not, the highly significant changes in profile 
shape above and below this value cannot be ignored. 



well modeled with a Gaussian function. It is only observed 
in Pop. B sources where it shows a redshift (FWHM in the 
range 1000-5000 km s~^). This can be the dominant com- 
ponent in some Pop. B sources. 

(iii) A BLUE component (blue shifted or blue asymmet- 
ric profile) seen in the H/3 profile of some Pop. A sources 
with i?Fcii= I(FellA4570)/I(H^) >0.5. It is more prominent 
in ClvA1549 where it can be defined as the residual emis- 
sion after subtraction of an unshifted Lorentzian BC in Pop. 
A sources. In this paper we model this component with a 
Gaussian. 

If there is a semi-coherent structure that gives rise to all 
or most of the broad line emission then all broad lines in a 
source should show similar properties allowing for some de- 
gree of stratification. H/3 apparently does not show the same 
compon ents in all sources mo ti vating the Pop. A-B nomen- 
clature jSulentic et al.ll2000al lbl. 120071 1 as an effective way to 
highlight this difference. We assume that the distribution of 
sources in the optical plane of 4DE1 has physical meaning in 
the sense that source occupation is driven by both physics 
and line-of-sight orientation. This assumption is reenforced 
by the discovery that sources in different regions of the op- 
tical plane (and especially Pop. A and B regions) show dif- 
ferences in almost all other multiwavelength measures. The 
sources chosen for this study allow us to test the hypoth- 
esis that the same three components observed in H/3 and 
Civ\1549 can account for the diverse line profiles of Lya, 
CivA1549, HeiiA1640, Sini]A1892, AhiiA1860, Cni]A1909, 
MgnA2800, H/3 and Ha. After testing this hypothesis ( ^4A[ 
we measure line intensities of the components ( i]4.3p . Anal- 
ysis of the properties of each line component yields insights 
about differences in physical conditions of the three identi- 
fied emitting regions. ( i]4.4|) . 



1.4 Defining an Empirical Approach 

This paper compares all of the principal broad lines (Lya, 
CivA1549, Cni]A1909, MgnA2800, H^) for a sample of 
sources where: (1) suitably high S/N spectra are available for 
all or most of the principal broad lines and (2) 4DE1 spectro- 
scopic diversity is sampled. Our approach is empirical and 
uses the H/3 profile fit for each chosen source to constrain 
the fits of all other broad lines in that source. High and low 
redshift Pop. B q uasars show rather symmetric Lya and 
CivA1549 profiles (|Barthel et al.lll990l : ISachev et all 120041 ) 
which is in sharp contrast to the appearance of H/3 in Pop. B 
sources where a red asymmetry is common. In Pop. A the sit- 
uation is reversed with blue asymmetric Lya and ClvA1549 
and symmetric Balmer line profiles. Comparisons of H/3 and 
ClvA1549 profiles suggest that three emission components 
are present in some or all sources. 

(i) A "classical" BC (FWHM = 600 -5000 km s"^). The 
H/3 profile is best modeled with a symmetric Lorentzian 
function for Pop. A sources (FWHM < 4000 km s"^) and a 
Gaussian function for Pop. B (FWHM > 4000 km s"^). This 
component may be weak or, even, absent in some Pop. B 
sources. The bulk o f the optical Fen emi ssion appears to 
follow the BC (e.g., IMarziani et~alll2003bl '). 

(ii) A VBC (FWHM~10000 km s"^). H/3 is reasonably 



2 SOURCE SELECTION AND OBSERVATIONS 

Sources were selected to be representative for the spec- 
tral diversity in the 4DE1 context. We selected six sources 
that include one quasar in each of the six most popu- 
lated 4DE1 sp ectral bins (Al, A2, A3, Bl, B1+ and B1++; 
ISulentic et al.l l2002). Figure [1] marks the location of these 
sources in the optical plane relative to the distribution o f 
the brightest 400-1- SDSS DR5 quasars (|Zamfir et al.ll2010l '). 
We remark that the validity of the analysis presented in this 
paper is not restricted to 6 objects only. Each object has 
been chosen as representative of each bin in the El plane 
of Fig. [T] These bin s are occupied by 90% of \ow-z quasars 
(|Zamfiret aLlbOld ). While high S/N and moderate resolu- 
tion spectral coverage from Lya to H/3 exists for very few 
quasars we were fortunate to find representative sources for 
these six bins. UV coverage from Lya to MgnA2800 comes 
from the HST archive for these low z sources save for 3C 110 
where MgnA2800 is covered in our optical spectrum. A list- 
ing of the available HST spectra can be obtained at the HST 
archive. A list with the datasets actually employed to build 
the UV spectra can be obtained from the authors|j Optical 



^ It is available at the URL: http:// web.oapd.inaf.it/~marziani/ 
uvlog.pdf. 
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spectra come from iMarziani et al.1 (|2003bl ): a log of obser- 
vation can be easily retrieved from their Table 2. UV and 
optical observations are generally not simultaneous, with a 
time mismatch of less than 3 years in all cases save Mark 
478; in this last case, the UV observations were collected 
more than 5 years after the optical ones. 

We make use of composite spectra for each of these bin s 
l|Sulentic et al.ll2002l : iBachev et aLll2004l : IZamfir et aLlbOloD 
computed fr om the bright SPSS sa mple or from our spectro- 
scopic atlas l|Marziani et aLlbOOSbl ). The Clll]A1909 blend is 
unavailable for 3C110 and the region of MgllA2800 has not 
been observed in 3C 249. Ha spectra are only available for 
three sources. Our composite spectra are used as surrogates 
for these gaps in our source data. 

Table [l] identifies the selected sources where: Col. 1 - 
lAU code name; Col. 2 - common name recognizable by 
NED); Col. 3 - adopted redshift; Col. 4: B magnitude; Col. 5 
- Galactic B-band absorption from NED ; Col. 6 - the spec- 
tral type following ISulentic et al.l (|2002r ) and Col. 7 - ad- 
ditional information. Two sources have been reverberation 



mapped (RM lPeterson et al.ll2004). Ga lactic absorption val- 
ues Ab are from ISchlegel et aU (|l998l ): if -Rv = 3.1, then 
Ay « 0.756Ab. This value has been used with the task 
deredden of IRAF. No internal extinction correction was 
applied to our sources. 



3 DATA ANALYSIS 

Data used in this paper are wavelength and flux calibrated. 
Spectra and spectral measures were shifted to rest frame 
wavelength and flux using the redshift listed in Tab. [T] after 
correction for Galactic extinction. 

The regions around Lya, ClvA1549 (including 
HellA1640), the blend at ^ 1900 A (including AI111AI86O, 
Silll]A1892, Clll]A1909), MgllA2800, H/? (including the 
FenA4570 blend), and Ha were analyzed using the specf it 
task of IRAF (Kriss 19941 ). We used an Fell template based 
on an ESQ spectrum of I Zw 1 to correct for Fell contamina- 
ti on. This templa t e was a lready employed for data analysis 



— .Marziani et al.l (l2003bl lcl) with improve ments described in 
IMarziani et al.l ([2009!). It is basically the Bor oson fc GreenI 



11993) template with two major modifications: the Fen 
emission underlying H/? was set according to a model 
simulation, due to the difficulty of disting uishing between 
H/? p roper and underlying Fe llopt features (jMarziani et al.l 
120091 ). In the red, our spectrum shows a significant feature 
at 6446 A In any case the Fe llopt correction is so small 
that it would not affect any of our results concerning the 
Ha line shape. In the UV, Fell emission around MgllA2800, 
Clll]A1909 and ClvA1549 is corrected a pplying a template 
obtained from a CLOUDY JFerland et al.l 119981 ) simulation, 
assuming ionization parameter 10~^'^^, electron density 
■|^gi2.25 (,jjj-3^ Felll emission, especially prominent around 
the I900A lin e blen d, was modeled using a template build 
following iVestergaard fc Wilkes (2001). There is consider- 
able uncertainty regarding the intensity of the Felll feature 
at I914A; this feature is heavily blended with Cin]A1909 
making it impossible a reliable estimate of the intensity 
of both lines in Pop. A sources (i.e., in any case the iron 
spectrum is strong). 

The three previously described line components were si- 



multaneously fit along with a power-law continuum and ap- 
propriate Fen (and Fein around 1900 A) template solution 
employing x^ minimization techniques. The relative inten- 
sity of the components was left free to vary although we im- 
posed the condition that velocity shifts and widths roughly 
match the ones measured in H/3 and in ClvA1549 (BLUE 
component only). 



4 RESULTS 

4.1 Multicomponent Fits 

Fig. [2] shows the resultant fits to the line profiles of Lya, 
CivA1549-HHeiiA1640, 1900A blend, MgiiA2800 and H^ for 
each source. Ha is not included since it does not add in- 
formation to the line profile interpretation based on H/3. 
Residuals are shown below each panel. Table [5] reports the 
total flux in the ClvA1549 line, and the normalized inten- 
sity ratios of BC, VBC and BLUE component; individual 
sources are discussed in Appendix [X] The following results 
are consistent with tren ds visible in median c o mposite spec- 
tra for the same bins (|Sulentic et al.l |2002| : iBachev et al.l 
l2004l : IZamfir et al.|[2O10l ). 

Broad Component The BC dominates LIL (Balmer lines. 
Fen, MgnA2800) emission in Pop. A sources, while it be- 
comes less prominent in Pop. B. The fraction of BC H/3 
flux decreases for Pop. B sources from bin Bl to Bl^ and 
from Bl"*" to Bl^"'". The 4DE1 sequence can be described 
as a sequence of decreasing H/3 BC prominence and may 
even involve BC demise. Evidence exists for AGN with a 
very weak or absent BC, some showing no emission lines 
(e.g. BL Lacs) or only narrow emission lines (type-2 AGN). 
Others more relevant to this study show lines with strong 
VBC emission that can be mistaken for a typical type-1 
AGN (|Sulenticet al.ll2000d . e.g. PG1416-129). We suspect 
that 3C 110 (also discussed in Appendix [A]) may be an ex- 
ample of this situation. Rather than the two component 
BC-I-VBC fits shown in Figure [2] and the correct interpre- 
tation for 3C110 may involve a single redshifted VBC com- 
ponent. FWHM ~ 10000 km s~^ would be unprecedented 
for the BC. 

In Pop. A sources I Zw I and Mark 478 Cni]A1909 is 
closely blended with Feni emission leaving open the possibil- 
ity that BC Cni]A1909 may be very weak or absent (for this 
reason no Cni]A1909 ratios are reported in Tab. [3]). Nearby 
AlniA1860 and Sini]A1892 are also very prominent in these 
two Pop. A sources. Conversely, bin Al source Mark 335 
shows stronger Cni]A1909 emission and AlniAlSGO intensity 
more similar to Pop. B sources. 

Blueshifted Component (BLUE) Figure [3] provides the 
justification for including a blueshifted component in the 
Figure [2] fits by showing the residual that arises from a two 
component fit to ClvA1549 (BC+VBC) for all six sources. 
Pop. A sources show a strong residual reflecting the well 
known blueshifted/asymmetric component. BLUE is present 
in the Lya and CivA1549 emission of all sources but 3C 110, 
with typical shifts and widths Ai; ~ —2000 km s~^and 
~7000 km s~^. Pop. A and B sources show similar shift 
and width values. Only bin A3 source I Zw 1 may show a 
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possible hint of this blue component in the H/3 line. BLUE is 
not seen in MgllA2800 while the 1900 A blend is too complex 
to allow any conclusion to be made. 



Very Broad Component A VBC component is one of 
the defining properties of the broad line spectrum in Pop. 
B sources. A VBC is not observed in NLSyl or even less 
extreme Pop. A quasars (see th e cautionary note re gard- 
ing very high luminosity quasars iMarziani et al.l 12003 ). Us- 
ing H/3 as a template we therefore fit only two components 
(BLUE + BC) in Pop. A sources. The residuals of those fits 
show no excess flux on the red wing of H/9 that might be 
interpreted as a weak VBC signature. This reenforces the 
earlier proposed ideas that FWHM~4000 km s~^ reflects 
some fundamental change in the structure and/or kinemat- 
ics of the line emitting region-perhaps connected to a critical 
value of the Eddington ratio. 

Restricting attention to our three Pop. B sources that 
show a VBC, shifts and widths are ~2000 km s~^ and ~ 
10000 respectively (3C 110 has two entries in Table [S] the 
second one refers to the assumption of no BC presence in the 
emission lines; see Appendixl^. These values are consistent 
with the ones derived for median spectra of more than 150 
Pop. B sources in the SDSS bright sample (Za mflr et al.l 
l2010f ). Using H/3 as the model we conclude that a VBC is 
also present in the ClvA1549 profiles of Fairall 9 and 3C 
249.1, giving rise to the extended red wings visible in Figure 
[3] These profiles look deceptively symmetric but both our 
empirical fits argue that ClvA1549 counterparts to the H/3 
VBC should be present. 



4.2 The HellA1640 profile 

The equivalent width of He""" lines is an Eigenvector 2 pa- 
rameter inv olving an apparent anti -correlation with source 
luminosity (JBoroson fc Greenlll993 '). The HellA1640 profile 
of Pop. B sources is customarily flat-topped and rather 
extended on the blue side. The proflle can be inter- 
preted as due to the sum of the same components ob- 
served in ClvA1549. An unambiguous decomposition of the 
HellA1640-|-Olll]A1663 blend for Pop. A sources is diflicult. 
A fit of the HellA1640-|- ClvA1549 profile reproduces the 
shape of the entire blend assuming that HellA1640 is due 
to blueshifted -I- broad component. In Pop. B sources, the 
analysis is even more complex because HellA1640 is blended 
with the red wing (VBC) of ClvA1549. The presence of 
blueshifted and VBC HenA1640 emission creates the ap- 
parent flat extension on the red side of ClvA1549 that 
has not been satisfactorily explained as yet. This exten- 
sion is visible in the sp ectra of several RL sources at both 
low- and high- redshift l|Marziani et al.lll996l : iBarthel et al.l 
1 199(3 ). This "plateau" effect is especially well seen on the 
red side of ClvA1549 in 3C 249.1 (Fig.g)) where the best fits 
for HenA1640 are obtained assuming the blueshifted and 
VBC components as for ClvA1549 and Lya. The same in- 
terpretation seems to apply to HenA4686: a good support- 
ing case involves the low- luminosity quasar PG 1138-1-22 
iMarziani fc Sulentid[l993 ). 



4.3 Intensity Ratios 

Tables [3l [4] [5] provide flux ratios for the three line compo- 
nents in each of the six sources. The rest-frame equivalent 
width for Lya is reported for each emission component in the 
last column of each table. Ratios are reported with 2 signifi- 
cant digits to avoid roundoff errors. Uncertainties on ratios 
are estimated from the errors in the individual component 
intensity which are typical ±20% for the strongest compo- 
nents (the ones of Lya and CivA1549, MgnA2800 BC, H^BC 
in Pop. A) or for individual or unblended features (Fell) and 
±50% for the weakest or most blended features (HenA1640, 
Sini]A1892, MgnA2800 VBC). Applying simple error propa- 
gation yields, for the intensity ratio of two components, la- 
beled 1 and 2: Aih/h) = ih/Ii)-^/{Ah/h)^ + {Ah/hY. 
For a typical error ~ 40% on line intensity, the error on ratio 

is A(/2//2)=0.56-(/2//l). 

BC The "classical" and almost ubiquitous broad line com- 
ponent shows some line ratios that yield important con- 
straints on the p hysical conditions: Lya/H/3 ~ 10 (see also 
iNetzer et al.lll995 ). compared to an expectation of ~35 for 
pure recombination; MgllA2800/Lya~ 0.1 - 0.3; relatively 
large 7?Feii(> 0.3). The Fell/H/3BC values assume that all or 
most Fell arises in the same region as H/3 BC. Removal of 
the VBC part of H/3 in the -Rpdi estimation increases val- 
ues over the previously published ones where no H/3 VBC 
correction was made, yielding -RfoII~ 0.5 for Fairall 9 and 
3C 249.1. ClvA1549/LyQ estimates range from ^ 0.15 for I 
Zw 1 and Mark 478 increasing to ~ 0.5 for Mark 335 and 
all Pop. B sources. An apparent anti-correlation between 
ClvA1549/Lya and -Rfoii as well as a positive correlation 
between Siiii]A1892/CivA1549 and Rfcii for Pop. A sources 
arc visible in median composite spectra of iBachev et al.l 
(2004). Our data in Table |3] confirm a positive trend be- 
tween Silll]A1892/ClvA1549 and i?Fcii for the 3 Pop. A 
sources, although low Silll]A1892/ClvA1549 values are ob- 
served for Efoii~0.5. The largest Siiii]A1892/CivA1549 
(and lowest ClvA1549/H/3) ratios are associated with the 
largest -RfoII Qt 0.5) estimates, probably implying metal en- 
richment ( ^4.4.1|l . 

BLUE The blueshifted component is less constrained due 
to its weakness in many lines/sources. The only ratios that 
can be measured involve Lya, ClvA1549, and, with greater 
difficulty HellA1640. Given the importance of the Lya/H/3 
ratio we made a special effort to derive an upper limit 
for it. They are estimated assuming that any H/3 emis- 
sion is peaking at less than 3cr the noise level, meaning 
that the upper limit to the line flux can be written as 



27ra « 1.0677p -FWHM, where 7p is the peak line inten- 
sity. The lower limits indicate that BLUE values of Lya/H/3 
are much higher than those of BC and VBC. This has im- 
portant physical consequences ( ^4.4[) . 



VBC Visually, the observed Fe llopt emission is fully con- 
sistent with an origin in the BC (i.e., no obvious evidence 
for any any Fen VBC). Even if we measure a FWHM of 
Fen consistent with the one of H/3bc for our 3 Pop. B 
sources, the VBC large width > 10000 km s~^ could cre- 
ate a pseudo-continuum underlying H/3. To settle the issue 
of a possible Fen VBC we attempted to fit a VBC of Fe 
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iiopt to t he Bl and Bl ^ com posites from the bright SDSS 
sample of IZamfir et~all (|201Q . S/N >200). The assumption 
of Fell VBC with the same shift and width of H/3vbc leads 
to implausible results with very large x^ unless Fen VBC is 
negligible. We conclude that we have no strong evidence of 
a significant Fen VBC emission. 

The red wing of the A1900 blend is ascribed to 
Clll]A1909 VBC emission, as the iron spectrum of Pop. B 
sources is always weak. 



i.3.1 The Mg\l\2800 /RI3 ratio 

With the goal of verifying whether the values of 
MgnA2800/H/3 reported in Table [3] and [5] are typical, we 
selected a sample of quasars brighter than g =19.0 with 
SDSS spectra that included both H/3 and MgnA2800 (red- 
shift range z = 0.40 - 0.75). While most of these spectra 
show moderate-low S/N they are good enough to allow one 
to assign the sources to 4DE1 bins in which context high S/N 
composite spectra can be computed (normalized to specific 
flux at 5100 A). We obtained 160 sources in bin Bl and 58 
in B1+. The VBC MgnA2800/H/3 ratios we derive from the 
bin Bl and Bl^ composite spectra are 0.7 and 0.5, respec- 
tively. Given Pop. A values in Table [S] and the consistent 
value MgnA2800/H/3~ 1.5 computed on the median spec- 
tra, we conclude that the VBC shows MgnA2800/H/3 ratios 
2-3 times lower than the BC. Due to the lesser prominence 
of the VBC in Pop. B sources, the FWHM o f MgiiA2800 
becomes significantly lower than that of H/? (jWang et al.l 
l2009l ). 



4.4 Physical Conditions 

We generated a multidimensional grid of CLOUDY 
(jFerland et al.l Il998l l simulations in order to infer ion- 
ization parameter U and electron density no values from 
our spectral measurements. Simulations span the density 
range 7.00 < log n^^ 14.00, and -4.50 =^ logt/ < 00.00 
with an interval o f 0.25 in both log U and logrio (c.f. 
iKorista et al.l 119971 ). Each simulation was computed for 
a fixed ionization parameter and density assuming plane 
parallel geometry. Part or all of this 2D grid of sinmlations 
was repeated assuming iVc= 10^\ 10^^ 10^^, 10^", 10^^ 
cm~^. Metallicity was assumed to be either solar or five 
times solar. Two alternative input continua were employed: 
(1) one assumed to be the standard AGN continuum by 
CLOUDY which is equivalen t to the continuum described in 
iMathews fc Ferlandl (119871) and (2) the low-z quasar contin- 
uum of iLaor et al.l (|l997al ). We exploit the simulations to 
deduce constraints on [/, A^c, and n^ from the most reliable 
intensity ratios we derived in this paper. 



BC The BC is the better constrained component with 
relatively large i?Fcii and ratios AhnA1860/ Sini]A1892, 
Siin]A1892/ClvA1549 requiring high density (log We ~ 12), 
low ionization (—3 < log (7 < —2) and large column den- 
sity (Ac> 10^"^cm~^). These conditions yield values of the 
important Lya/H/3 and MgnA2800/LyQf ratios that are con- 
sistent with the ones observed, and explain -Rfcii < 0.5 — 1.0 
as well. Larger -Rfcii values suggest super-solar metallicity 
fi ]4.4.1|) . Hereafter we will refer to the BC-emitting region as 



the LIL-BLR while we will retain the term BLR in its more 
general meaning i.e., as the region where all of the broad 
emission is produced. 

BLUE The very large Lya/H/3 ratio suggests radically dif- 
ferent conditions including lower optical depth and consis- 
tency with 7Vc< 10^-^ cm"^ (see '^^ . The Lya/H^ and the 
HenA4686/ ClvA1549 ratios suggest high ionization log U > 



-1 and densities in the range nc~ 10 



-10^ 



cm ^. There 



is no detection of a Cm] A1909 blueshifted component in the 
A1900 blend, but this is consistent with high ionization gas 
at density logne~ 10: the expected Cin]A1909/ClvA1549 ra- 
tio for logU ~ —1 is ^0.1, and any blueshifted Clll]A1909 
emission would be lost in the 1900 A blend. We included 
a component with the same shift and width derived for 
ClvA1549 in the fits of F9 and 3C 249.1, but specf it yielded 
or negative values in the final fit. 

VBC If i?Feii~0 in the VBC and given the low 
MgnA2800/LyQ ratio, we infer a high ionization parameter 
log [/ ^ — 1. However, Results for VBC are contradictory 
with all ratios consistent with high ionization but some ra- 
tios indicating high and others low density. 

We observe Cni]A1909/ClvA1549~ 0.1 - 0.2 which in- 
dicates a relatively low density emitting region with nc~ 
lO'J - 10^0 cm-^ The observed HenA4686/ClvA1549 ra- 
tio (0.1 - 0.2) is also consistent with this picture. Much 
higher density would collisionally quench Cni]A1909 to un- 
detectable levels. Not all VBC line ratios are consistent with 
moderate density picture. If the Balmer decrement is very 
flat, and the Lya/H/3 ratio in between 5 and 10, then very 
high density gas (ne~ lO^^cm"^) is required (jNetzer et al.l 
[1995) as well as a very large optical depth if photoionization 
is assumed. A large column density Ac> lO^^cm"^ yields 
condi tions consistent with these line ratios measured for the 
VBC. ISnedden fc Gaskelll (|2007l ') computed line profile ratios 
Lya/Ha and Ha/H/3 for several luminous Seyfert 1 nuclei 
(mainly Pop. B) and reached similar conclusions in favor of 
large We and A'c. 



4-4-^ Chemical Abundances and Fein Contribution 

The LIL-BLR is the only region allowing for tentative abun- 
dance considerations; results should obviously be applicable 
to the other emitting regions. Ratios involving metal lines 
should be considered the most robust for deriving physi- 
cal inferences. Intermediate-ionization lines (like AlinA1860 
and Silll]A1892) are produced in a fully ionized region 
which ma kes them less depend ent on column density than 
the LILs (jNegrete et al.M2010h . From the simulations, we 
find that the ratio AhnA1860/ Siin]A1892 is not strongly 
dependent on Z: the ratio increases by about 40% pass- 
ing from Z — IZq to Z — 5Zq, for lognc~ 12 and 
logU «-2. However, CivA1549/ Siin]A1892 (or, equiva- 
lently ClvA1549/AhnA1860) is much more affected. The low 
Sini]A1892/ClvA1549 and the other line ratios we measure 
for Mark 478 can be all accounted for with Z = 5Zq, log 
nc~ 11.6, and \ogU ~ -2.4. No solution with a well defined 
value of We and U is found ii Z — Zq is assumed. Even higher 
metallicity, ~ WZq, is strongly suggested for I Zw 1. The 
assumption of metallicity above s olar is consistent with ob- 
servations of high-z quasars (e.g. iHamann fc Ferlandlll999l : 
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I Juarez et al.ll2009h . and the two sources where the evidence 
of higher abundance is strong also show the most prominent 
NvA1240 hne in the Lya panels of Fig. [21 as expected. 

In the cases of high metallicity we expect that Feiii 
should be especially prominent in the 11900 spectral region. 
The physical conditions derived for the LIL-BLR remind 
to the Weigelt blobs of r] Carinae, located in the equatorial 
plane of the system, perpendicular to the symmetry axis of 
the bipolar lobes forming the "Homunculus" nebula. The 
Weigelt blobs are believed to be dense gas photoionized by 
the radiation associated to the central, massive star and to a 
possi ble companion (e.g., [Johansson et al.ll2000l : iDavidsonI 
120051 ). The spectrum of the Weigelt blobs shows very weak 
ClIl]A1909 along with a prominent line at 11914, ascribed to 
the z'^P^ — ^ a^Ss Felll transition. Although this is not a res- 
onant line (the lower level is « 3eV above ground), the line 
appears very strong because the upper level is populated by 
Lya fluorescence. Indeed, in I Zw 1 where lines are narrow, 
the peak emission at 11914 is actually visible, and the Felll 
e mission around 1900 A may r esemble the one of Model 3 
of lVestergaard fc Wilkei l|200iri . 



5 DISCUSSION 

In the following we will discuss measurements in fixed in- 
terpercentile velocity ranges ([J5TT} as an approach comple- 
mentary to our multicomponent analysis. After consider- 
ing properties of the BLUE component ( ij5.2|) we will show 
how our results on the physical properties of the different 
emitting regions do not favor a predominance of line emis- 
sion from an accretion disk f H5.3.ip or from a binary BLR 
f q5.3.2|) . We suggest that the properties of the BLUE, BC 
and VBC emitting regions as well as their trends along the 
El sequence can be explained by the interplay of gravita- 
tional and radiation forces ( ^5A\ and i]5.5p . Last, we will 
discuss our interpretation of the BLR and compare it to 
past work f> i5.7p . 



5.1 Interpercentile Velocity Measurements 

The true profile shapes of the various line components are 
uncertain. The Lorentz-like shape of Pop. A II/3 profiles can 
be considered the most accurate assuming that the peak 
and the wings do not represent ind ependent emission com- 
ponents (e.g., iLa Mura et al.ll2009l ). The lack of inflections 
in ev en the highest S/N composite profiles (IZamfir et al.l 
l20ld ) motivated us to favor the single component Lorentz 
fits. Uncertainties in our decompositions force us to restrict 
the analysis to estimates of the relative component inten- 
sities, average line shifts and widths. A double Gaussian 
function is a good approximation but the individual Gaus- 
sian fits may provide uncertain estimates of the relative BC 
and VBC fluxes in the core of the line where the two compo- 
nents overlap in radial velocity. It is difficult to say whether 
a stronger core component and a less broad, fully redshifted 
VBC provide a wor s e fit, a lthough the SDSS composite spec- 
tra oflzamfir et al.l (|2010| '1 and earlier work on PC 1416-129 
l|SuIentic et al.l l2000d ) favor our BC-VBC decompositions. 
Also, if our interpretation of some sources showing II/3 that 
is almost pure VBC H/? (3C110 in this paper) is correct. 



then the Gaussian functions that we employ represent a rea- 
sonable description of the VBC profile shape. Measuring line 
ratios at fixed inter-percentile velocities is an alternative ap- 
proach that might have a more direct physical meaning if the 
BLUE and VBC matter only in the highest radial velocity 
ranges. We therefore report emission line fiuxes measured in 
three radial velocity ranges of Av^— 2000km s~^ centered 
at km s"^ ("0" range), ±4000 km s"^ (B and R ranges 
for Pop. A sources) and km s~^, ±6000 km s~^ (Pop. B 
sources) in Table (6] The reported measurements were car- 
ried out after subtraction of Fen, Fein and of contaminating 
narrow and broad lines. Errors are estimated at a 2a con- 
fidence level as twice the rms scatter in the radial velocity 
range times the AA in A. Two different centerings of the ra- 
dial velocity ranges are needed to account for the line width 
difference of Pop. A and B sources. The overall interpreta- 
tion as deduced from the multicomponent fit is confirmed. 
In particular, the data of Table [B] indicate that there is no 
evidence of an additional red component for Pop. A objects. 
There is also no convincing evidence of a BLUE component 
in II/3 of Pop. B sources. The large error bars associated 
with the Lya/H/S ratios support the same conclusion moti- 
vated by the multicomponent fits for the BLUE component: 
the upper limits and the possible detection of I Zw 1 imply 
LyQ:/H/?3> 10. Line ratios for the B, and R ranges yield 
the same qualitative conclusions as the multicomponent fits, 
with the limitation that we cannot obtain meaningful results 
for the heavily blended features at A1900. Also, the large er- 
ror bars in the H/? R range of Pop. B sources refiect the 
residuals of [O ill] 14959 subt raction. A compar ison of our 
Lya/H/3 ratios with those of (jNetzer et al.lll995l ') and of our 
H q/H/3 ratios (from t he m ulticomponent fits) with those 
of ISnedden fc Gaskelll ((2007]) show agreement especially in 
the redshifted line wings i.e., for the VBC. This suggests 
that our approach yields consistent values although the true 
shape of the VBC component is uncertain especially in the 
range -1000km s"^^Ur^ 1000km s"\ 

5.2 Optically Thin Gas? 

The weakness of the response to continuum changes in 
the H/3 wings and inter-line comparisons e.g. have mo- 
tivated arguments th at the VBC is mainly produced 
in optically thin gas JMorris fc Ward! Il989l: IZhend Il992l: 
Shields et al] 1 19951 : ISulentic et al.l l2000d : ICorbin fc SmithI 
20001 ). On the other hand, several authors have inter- 



preted the response of the blue wings in ClvA1549 and 
H/3 as evidence for predominan t ly optically thick gas 
(e.g. iTurler fc CourvoisieJ Il997l: iKorista fc GoadI |2004| : 



ISnedden fc Ga skcU 2007; Shapovalova et al.ll20ld . and ref- 
erences therein). Korista & Goad (2004J) showed that an op- 
tically thick VBLR can explain the stronger core response of 
H/3 to continuum changes. The small Lya/H/? and Ha/H/3 
ratios suggest a large column density A^c for both LIL-BLR 
and VBLR. This interpretation may be correct for some 
transient sources but is unlikely to apply to more stable 
sources that show a strong red-ward asymmetry (i.e. Pop. 
B quasars). In an optically thick medium the intensity of a 
recombination line depends on the luminosity of the ionizing 
continuum; if the medium is optically thin it depends on the 
amount of emitting gas. In practice this places a constraint 
on the column density A'^c. In high luminosity quasars op- 
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tically thin gas may not be adequate, even with a covering 
factor approaching unity, to explain the equivalent width 
and luminosity of the ob served VBC (JMarziani et alj 120061 : 
ISnedden fc GaskellllJOOTl ). 

Is the blueshifted component due to optically thin gas 
perhaps ejected from the central engine? Fig. [S] shows the 
predicted equivalent width of Lya and Lya/H/? intensity 
ratio (thin lines) as a function of the ionization parameter 
U for iVc= 10^^ cm~^ and iVc= 10^^ cm"^. Computations 
Eissume total covering, solar metallicity and ne= 10^°cm~''. 
W(Lya) remains approximately constant until the optically 
thin domain is entered at a "critical" value of U . Beyond that 
value W(LyQ) is roughly proportional to the inverse of U: in 
the optically thin domain the increase in continuum lumi- 
nosity does not give rise to a corresponding increase in line 
luminosity. The most stringent observational boundary is set 
for Mark 478 and Fairall 9, with LyQ/H/3>45. At the same 
time, the median CivA1549/LyQ and HenA1640/CivA1549 
values suggest —1.5 < f/ ~ —1 if logne= 10. Considering 
that the observed W(Lya) is in the range 15 - 60 A, low col- 
umn density gas at A'c= 10^^ cm~'^ can explain the largest 
observed equivalent widths only if there is substantial cov- 
ering of the source, i.e., the source is almost completely cov- 
ered by a geometrically thin shell (~ 10^^ cm) of optically 
thin gas. A smaller covering fraction is required if A'c= 10^^ 
cm~^ or larger. 

5.3 Alternative Interpretations 

5.3.1 Accretion Disk Contribution 

Some Pop. B sources show very broad double-peaked H/3 
profiles that have been interpre ted as the signature of ac- 
cretion disk line em ission (e.g., lEracleous fc Halpernll2003l : 
IStrateva et al.luOOa 'l. A double peak ed signature has bee n 
found for MgiiA2800 in Arp 102B (|Halpern et all 1 19961 ). 
Others have argued that double-peaked profiles are too 
rare and disk model soluti ons too disparate to justify the 
disk emission hypothesis (|Sulentid 1 1989). Attempts have 
been made to find lar ger samples of double-peaked sources 
IStrateva et a l. 2003; Bon et al.ll2009l '): however, most (90%) 
of these source are not double-peaked but rather red or 
blue asymmetric. A large sample of t he highest S/N SPSS 
spectra binned in the 4DE1 context (jMarziani et al.ll2009l : 
IZamfir et al.ll201Cl ) do not show any hint of a double-peak 
signature in high S/N (~250) composites. Broader pop. B 
sources unsurprisingly favor red asymmetries. Sources in- 
cluded in a given Pop. B bin of Figure 1 show a diversity 
of red/blue bumps and asymmetries however bin composite 
spectra show a BC-I-VBC blend similar to the one observed 
in e.g. Fairall 9 and 3C 249.1 of this paper. Where is the 
disk signature? If a significant part of the broad lines is 
produced in an accretion disk then it must be obscured or 
mask ed by non-disk emission components (e.g. IZhang et al.l 
[2009'). For Pop. B sources, one could argue that the sum 
of the BLUE and VBC represents a double peaked com- 
ponent accounting for the broadest part of the line profile 
with the BC emitted elsewhere. A difficulty with this inter- 
pretation is that the gas in the BLUE and VBC appear to 
be arise in regions with different physical conditions. The 
BLUE component is not consistent with high density and 
large column density as expected in standard accretion disk 



models (e.g., IColIin-SoufFrin fc Dumontlll990l l. In the case 
of Pop. A sources the narrower Lorentz-hke profiles, inter- 
preted as a disk signature, would require line emission from 
a disk with radius ~10^ gravitational radii which is theoret - 
ically and observationally disfavored (| Collin fc Zahnlll999l ). 
The evidence continues to grow weaker for a predominance 
of accretion disk emission in Pop. A and Pop. B line pro- 
files. Only a minority of very-broad sources whose proto- 
type is A rp 102B remain consi s tent with accretion disk pre- 
dictions JStrateva et all l2003l : lEracleous fc HalpernI l2003l : 
ICezari et al.ll200'7^ . 



5.3.2 Binary BLR 

If quasar activity is driven by merger and accretion pro- 
cesses then we might expect to find some/many sources 
with evidence for a binary black hole. At least some 
of these sources might involve two AGN with associ- 
ated BLRs. It has been proposed that binary BH con- 
figurations mi g ht give rise to double - peaked line profiles 
JCaskelJ Il985l: iHalpern fc Filippenkol Il988l: IShen fc Loebl 



I2OO9I : iLauer fc Borosonll2009l . see lGaskellll2009al for a recent 
review). Two co-orbiting BH each with a BLR producing 
(typical) single-peaked broad lines could give rise to double 
or single peaked composite line profiles depending on the 
line-of-sight orientation and radial separation of the BHs. A 
large peak shift v^ > 2000 km s~^ is not unphysical since 
Vt ~ 2000 km s~^ is less than the typical virial velocity at 
the reverberation estimated radius of the BLR. 

This idea is not easily applicable to Pop. A sources 
which show quite narrow and symmetric H/3 profiles. The 
situation for Pop. B sources appears more favorable if 
one considers, for example, a binary BLR origin of the 
BLUE and redshifted VBC components. The reductio ad 
extremum would view all pop B sources as binary pop A 
sources. One problem with this view would involve the fact 
that virtually all pop. B sources are pop. B. However in 
some pop. B sources unshifted BC gas might be associ- 
ated with a low-ionization region at larger radius and dis- 
tributed around the center of mass of the two orbiting 
BH. The radial velocity of the first component would be 
^r,i ~ ^/GMjd ■ sini ■ sin(27ri/P) • m^jM^ while for the sec- 
ond i;r,2 = ^GM/d ■ sini ■ sin(27rt/P -I- tt) • mi/M at any 
given orbital time t. P is the orbital period, ttii and 7712 the 
BH masses, with M = mi +m2, d the black hole separation. 
The inferred properties of the BLUE and VBC components 
lead one to infer that the masses of the two BHs should 
be roughly equal. The Mbh similarity again raises questions 
about why we find very different Lya/H/? ratios and inferred 
physical conditions for the BLRs around the two BH. In ad- 
dition, the broadest profiles should "narrow" dramatically 
on a time scale of « 7rd/2ii ~ lOd iefJnnnyr- Note that the 
profile presented in Figure 10 of the lShen fc Loebl ([20091) pa- 
per is not double-peaked but a prototypical Pop. B sources 
showing the red asymmetry that we suggest is due to a dis- 
tinct VBC emitting region. If one ascribes a (broader) core 
component to a more massive black hole and the redward 
asymmetry to a narrower component (less massive BH), then 
the signs of the two component radial velocities Ur,i and iir,2 
should always be opposite, which is not the case in the ma- 
jority of Pop. B sources. The expected change in radial 
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velocit y is not seen in Arp 102B-like double pea ked sources 
either JEracleous et al.lll997l : [Gezari et al.ll200i1) . 



5.4 Role of Mass and Eddington Ratio 

Table [7| reports the Mbh and -L boi/I/Edd values for the 
source s following the prescription of lVestergaard fc PetersonI 
l|2006r ) where: Col. 1 - source name; Col. 2 - specific fiux at 
5100 A; Col. 3 - bolonietric luminosity assuming a correc- 
tion factor of 10 from the observed luminosity at 5100 A; 
Cols. 4-5 - FWHM of the entire H/3 profile and FWHM of 
H/3bc as derived from the x^ fits; Cols. 6-7 and 8-9 present 
the corresponding black hole masses and Eddington ratios. 
Mbh values are reduced by A log Mbh~ 0.3 when the BC 
component is employed. The effect is larger (factor of 10) if 
ClvA1549 is used as a virial proxy for H/3 in Pop. A sources 
ISulentic et al.l2007^ . Sources follow a sequence of increasing 
mass and especially decreasing I/boi/^Edd. Pop. A sources 
show the highest Eddington ratios by an order of magnitude 
or more. A correction by radiation pressure would leave the 
Eddington ratio trend unaffected (jMarconi et al.ll2008l V It 
is interesting to note that the VBC occurs only in Pop. B 
sources where I/boi/iEdd< 0.1, but BLUE is strong in both 
Mark 478 and 3C 249.1 whose I/boi/iEdd is different by al- 
most an order of magnitude. 



5.5 Gravitational and Radiative Acceleration 

If line emitting gas is optically thick to the Lyman contin- 
uum then the radiation will exert an outward acceleration 
that is inversely proportional to the column density and pro- 
portional to the ionizing luminosity. The ratio between the 
radiative and gravitational accelerations is: 



flrad 



-1 

^BH,8-' 



(1) 



where Mbh is in units of 10 solar masses, and L44 is the 
luminosity of the ionizing continuum (A < 912A) in units of 
10*'*ergs s"'^. The equation can be written in the convenient 
form 



Q^rad 



7.2 



Lhoi 



n: 



(2) 



If ''a ^ 1 then radiative acceleration dominates. This may 
be the case for the BLUE component if we interpret the 
blueshift as a Doppler shift. A condition of equilibrium may 
be reached at ra ~ 1. If I/boi/iEdd= 1, the corresponding A'c 
is ^ 10^* cm~'^, a plausible value for the LIL-BLR. Finally 
if ''a «C 1 the emitting gas may be unable to withstand the 
central black hole gravity and may fall toward the center 
giving rise to the observed redshifted VBC. This interpreta- 
tion is compatible with large A'c values for the VBLR and 
would naturally explain why the VBC is observed in objects 
with low I/boi/^Edd ~ 0.1. Of particular relevance to Pop. 
A sources are emitting clouds moving under the combined 
effect of gravitation and radia tive acceleration (jMathewd 
1 19931 : iNetzer fc Marzianil l20ld ) that produce the observed 
Lorentz-like profiles. If I/boi/iEddJ;0.1, as in Pop. B, gas of 
log A^c~ 22 would be pushed away while gas even larger 
column density could flow out in Pop. A sources. 



5.6 Radio Loudness 

If we compare the H/3 profile of RQ and RL in Pop. B 
i.e., where most radio loud objects are found, there is no 
major difference although RL objects show a preference to- 
ward more extreme red ward asymmetries (jMarziani et al.l 
[2003g). If we look at high ionization lines of RL objects, the 
C1VAI549 profile can show a blueward asymmetry although 
a redward asymmetry is more frequ ent (Marziani et al. 199y; 
ISulentic et all l2007l : IPunslvl I2OI0I '). Radio quiet objects of 
Pop. B show a somewhat different behavior: ClvA1549 
BLUE is systematically stronger and the pro file is more of- 
ten symmetric or blueward asymmetric. The lSulentic et al.l 
(2007 ) data suggest that radio loudness may modulate the 
relative intensity of BLUE and VBC, lowering BLUE, per- 
haps because radio ejections aff ect the radial flow seen in RQ 
quasars (.Marziani et al.l 119961). BLUE is however deflnitely 
observed in RL objec ts (|Barthel et al.l Il990l : IWills et al] 
ll995l : ISmaUet al.lll997h . 



5.7 Previous Studies and the Present 
Interpretation of the BLR 

We are aware of only one other attempt to m odel all ma- 
j or lin e proflles in a sample of low z quasars. IShang et al.l 
(|2007l) studied 22 sources and modeled all lines using two 
components-broad and very broad. Their sample is larger 
and allows for some correlation analysis. However, Pop. A 
sources are likely overrepresented, a result of their use of 
an X-ray select i on th at favors Pop. A sources. The results 
of IShang et al.l ((20071) should be therefore related mainly 
to three Pop. A sources included in our investigation even if 
some Pop. B present in their data suggest a change in trends 
at FWHM(H/3)Ri 4000 km s"^ (their Fig. 5). Their major 
empirical result, namely that ClvA1549 and Lya are very 
different from Ha and H/3 in terms of line width, asymme- 
try, and velocity shifts is confirmed by the present study as 
far as H/3 is concer ned. Our results also suggest, along with 
IShang et al.l ((20071), that a simple radially stratified ioniza- 
tion structure of the BLR is unable to explain the different 
line profiles. A wind component has been isolated through 
our El-based approach and three components appear to be 
a minimum set for the reasons summarized in the introduc- 
tion. The 4DE1 approach unifies the lines within specific 
sources and allows for the definition of a trend among the 
sources. BLUE and LIL-BLR are characterized by well de- 
fined properties. The VBLR is a region that could be defined 
by virtue of high ionization and, if our considerations are 
correct, large A'c. Within these constraints a large range of 
parameters - especially density - is suggested. The VBLR is 
therefore the only region where a locally optimally emitting 
cloud sch eme - where a range of properties occurs at every 
location ((Baldwin et al.lll995l )- could apply. 

Recent reverberation mapping studies attempted to 
consider resp onse of lines to contin uum chang es in bins of ra- 
dial velocity ((Bentz et al. 2008; D ennev et al . 2009). Infall is 
suggested by the observation of a faster response to contin- 
uum change in the red wing of H/3 in NGC 3516 and Arp 151. 
Arp 151 belongs to a perturbed system, and infalling gas can 
give rise to a redward asymmetries on spatial scales much 
large r than the one of the B LR (e.g., Rafanelli fc Marzianj 
[ 19921 : [Marziani et al.[[2003al ). Even if it is too early to gen- 
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eralize, these results are consistent with the interpretation 
of the H/3 redward a symmetry /redshift in Pop. B sources 
(e.g., iGaskeHlboOObl ). As far as the VBC is concerned, a 
possible problem with an infall interpretation is related to 
the line anisotropy. Very high A'c clouds should prevalently 
emit lines facing the centr al photoionizing con tinuum source; 
this is true also for H/3 JFerland et allliggd l. Therefore, if 
line emission is anisotropic, gas should be visible from the 
far side of the continuum source and infalling gas should 
yield a net blueshift. However, if emission comes from the 
non- illuminated side of the cloud (jFerland et al.ll2009l '). and 
the more distant part of the infalling gas is obscured a net 
shift to the red is expected as frequently observed in the line 
profiles. 



The low-ionization BC seems to be preferentially as- 
sociated with young, or rejuvenated, systems like NLSyl 
galaxies l|Mathuil I2OO0I : ISulentic et al] l2000al . |2008| ) whose 
iboi/iEdd is relatively large (0.2 <Lboi/iEdd< 1). En- 
hanced star formation in NLSyls ha s been recently d emon- 
strated from analysis of Spitzer data (ISani et al.ll 2010) while 
clues have existed for a long time (e.g. ISanders et al. 19881 : 
iKrongold et al.ll200ll . and references therein). Higher metal- 
licity is likely a consequence of top -heavy starb ursts that 
may be frequent in galaxy nuclei fe.g. iBonnell fc R ice 2008) ■ 
Such starbursts are expected to sustain the high accretion 
rate typical of Pop. A sources and to provide the enriched 
LIL-BLR material in extreme Pop. A sources like I Zw 1. 



5. 7. 1 On the Origin of the Broad Component 

The BC may be distinct from the other emitting regions 
and in many sources (Pop. A) may be the dominant or only 
visible component in all LILs. Two defining BC properties 
involve low ionization degree and symmetry of the line pro- 
files. We find BC width FWHM^ 5000 km s'^ith frac- 
tional peak shift Aur/FWHM < 0.1. This leads to the con- 
clusion that the LIL-BC may be the only H/3 component 
that is a valid virial estimator. The width of Fen and espe- 
cially of MgllA2800 likely offer the next best virial estimators 
l|Sulentic et al.ll2006l : I Wang et al.ll2009D . 

The necessity of a low-ionization emitting sub-region 
within the BLR has been s tressed since long, especia lly by 
S. Collin and collaborators l|Collin-Souffrin et al.lfl988l '). Ob- 
servational support in favor of this idea came by the real- 
ization of the very different profiles of LILs and HILs, as 
reviewed in the introduction. More recently, the analysis of 
the Call IR triplet suggests that a rather dense gas with 
density log rie 11.5 and log U ~ —2.5 is present in the BLR 
and that is primarily responsible for the observed Ol, Ca 
II, and Fell lines, based on the resemblance of their profiles 
(jMatsuoka et al.ll2007, 2008 ). There is a close agree ment be- 
tween the values derived bv lMatsuoka et aU l|2007l 'l and the 
values derived in this paper through t he ratios AI111AI86O/ 
Siiii]A1892, C1VAI549/ Siiii]A1892 (see lNegrete et al.ll20ld . 
for a more detailed report). 

LIL emission from AGN has been the subject 
of a controversy re l ated to t he ionization mechanism: 
ICoUiii-Souffrin et al.l (|l988l ) and ICoUin-Souffrin fc DumontI 
1 19891 ) suggested mechanical heating to explain the LIL 
prominence and especially the Fell intensity. Explaining 
strong LIL emission under the assumption of solar abun- 
dance has been a challenge for photoionization models. The 
ClvA1549/HellA1640 and Silll]A1892/ClvA1549 intensity 
ratios lead us to conclude that the assumption of metal en- 
richment is appropriate. Among Pop. A objects there could 
be an important role of metallicity in farther defining the 
properties of the LIL-BLR. We see indeed a trend in the 
prominence of the NvA1240: the most prominent NvA1240 is 
observed in I Zw 1 and Mark 478. If the NvA1240/ClvA1549 
ratio can be used as a metallicity indicator, then our six ob- 
jects might foll ow the accretion rate-metallicity relationship 
pointed out bv lShemmer et al.l ([2004). An increase in iron 
abundance from a few to sever al times solar (i n agreement 
with much previous work, e.g. iHamann fc Ferland 1999 ') is 
suggested for sources with 7?Feii> 1 (A2-|- bins). 



6 CONCLUSION 

Ignoring the composite nature of broad emission lines in 
quasars will likely result in a poor or erroneous understand- 
ing of their emitti ng regions (as known for 30 years; e.g., 
[Gordon et al.lll98lh . The contextualization of 4DE1 not only 
shows that quasars are spectroscopically diverse but it pro- 
vides a way to organize the differences in order to facilitate 
interpretation of the emitting components in different classes 
of sources. Analysis of 6 low-z quasars that explores 4DE1 
diversity and are believed to be representative of the wide 
majority of low-z quasars suggests that the strongest opti- 
cal/UV emission lines involve one or more of three different 
components that can be defined by their shifts with respect 
to the rest frame of each quasars: a redshifted VBC, an un- 
shifted BC component, and a blueshifted (BLUE) compo- 
nent. The three components have been tentatively isolated 
through models of the strongest lines from Lya to MgllA2800 
and Fe llopt using H/3 as a template, and characterized as 
follows. 

(i) The BC is the broad line component present in the 
overwhelming majority of type-1 AGN. Fell and MgllA2800 
emission likely arise in the same region. Other broad lines 
usually show a BC component although there may be ex- 
ceptions: Clll]A1909 may be very weak in extreme Pop. A 
sources where metallicity enhancement is most in evidence. 
At the other extreme of the 4DE1 sequence some Pop. B 
sources may show only weak or even zero BC emission. 

(ii) Pop. A sources show dominant BC emission in LILs 
and strong BLUE component emission in HILs. The BLUE 
component is most prominent in NLSyls, but is detected 
in most quasars, including Pop. B sources. The prominence 
of the BLUE component in extreme Pop. A NLSyls (al- 
most always radio-quiet) suggests that it is not closely re- 
lated to any relativistic outflow due to radio jets. It has 
been interpreted as the signature of an accretion disk wind 
where the receding p art of the flow is obscured by the disk 
(jMarziani et al.lll996r ). If the interpretation of i]5.5l is correct, 
then BLUE may be due to gas upon which radiative forces 
dominate. We are able to constrain ionization parameter, 
density and column density of this outflowing region. These 
quantities are important to assess feedback effects on the 
interstellar medium of the quasar host galaxy. 

(iii) The redshifted VBC component can be described as 
a deflning LIL property in Pop. B sources. The BC is dis- 
tinguished empirically from the VBC by its much larger 
FWHM and signiflcant redshift. The low Balmer decrement 
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and LyQf/H/3 ratio suggest, if we assume photoionization 
as the heating source of the gas, high ionization and large 
cohimn density. The physical conditions in the VBLR are 
not clear, and probably involve a stratification of properties 
that we are unable to resolve at this time. On the one hand, 
the Cin]A1909/ClvA1549 ratio suggests high ionization and 
moderate density. On the other hand, the strong HellA1640 
emission, the low LyQ/H/3 ratio as well as the flat Balmer 
decrement (Ha/H/3« 1) are better explained by extreme 
density at very high ionizatio n. In this region a loc ally opti- 
mally emitting cloud scheme IjBaldwin et al.ll 1993 1 could be 
appropriated. 
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Table 1. Selected Sources 



lAU Code 



NED Name 



"IB 



Ab Sp.T. Comments'^ 



J00535+1241 


UGC 00545 


0.0605 


14.4 


0.28 


A3 


= I Zw 1, Ha 


J14421+3526 


MRK 0478 


0.0771 


15.0 


0.06 


A2 




J00063+2012 


MRK 0335 


0.0258 


14.0 


0.15 


Al 


RM, Ha 


J01237-5848 


Fairall 9 


0.0462 


13.5 


0.12 


Bl 


RM. Ha 


J11042+7658 


PG 1100+772 


0.3116 


15.7 


0.15 


B1+ 


= 3C 249.1, no Mgll 


J04172-0553 


3C 110 


0.7744 


15.8 


0.19 


B1++ 


no CHI 



RM: Reverberation Mapping data; listed in lPeterson et al.l l|2004h . 



Table 2. Relative Intensity of BLUE, BC, and VBC in the ClvA1549 Line 



Sp. T. 


Name 


F(CivA1549;[^ 


BC(Civ)/ 


VBC(Civ)/ 


BLUE(Civ)/ 








F{Civ) 


F{Civ) 


F{Civ) 






[ergs s~'-'-cm^^] 








A3 


I Zw 1 


0.74E-12 


0.36 


0.00 


0.74 


A2 


Mark 478 


1.17E-12 


0.61 


0.00 


0.39 


Al 


Mark 335 


6.31E-12 


0.87 


0.00 


0.13 


Bl 


F9 


5.84E-12 


0.43 


0.27 


0.30 


B1+ 


3C 249.1 


2.11E-12 


0.17 


0.31 


0.38 


B1++ 


3C 110 


1.60E-12 


0.00 


1.00 


0.00 



F{Civ) = BLUE + BC + VBC 



Table 3. Broad Component (BC): 1000 km s-^< FWHM < 5000 km s"!, Av^ ~ km s" 



Sp. T. 


Name 








Intensity Ratio 








W^ 




Civ/ 


Aim/ 


Sim/ 


Mgii/ 


Fen/ 


Ha/ 


Lya/ 


Mgii/ 


Lya 






Lya 


Sim 


Civ 


Lya 


H/3 


H/3 


H/3 


H/3 




A3 


I Zw 1 


0.14 


0.61 


0.97 


0.31 


1.80 


4.2 


5.5 


1.70 


61 


A2 


Mark 478 


0.14 


0.41 


0.47 


0.09 


0.95 




10.3 


0.94 


87 


Al 


Mark 335 


0.45 


0.33 


0.07 


0.07 


0.47 


4.6 


15.8 


1.03 


102 


Bl 


Fairall 9 


0.53 


0.32 


0.18 


0.21 


1.04 


<6.1 


12.3 


2.54 


158 


B1+ 


3C 249.1 


0.44 


0.68 


0.17 




0.54 




10.7 




32 


B1++ 


3C 110 


0.55 






0.08 


< 0.36^ 




17.0 


1.3 


71 



" Equivalent width in A. 

Assig ned following the criterion that related FWHM — S/N and minimum detectable intensity, as in lMarziani et al.l 
ll2003bli 
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Table 4. Bluoshiftcd component (BLUE): FWHM ~ 7000 km s-\ Avr 



-3000 km s- 



Sp. T. 


Name 




Intensity Ratio 




w« 




Civ/ 


Hell/ 


Ha/ 


Lya/ 


Lya 






hya 


Civ 


H/3 


H/^i 




A3 


I Zw 1 


0.25 


0.41 


4.2 


~18 


60 


A2 


Mark 478 


0.66 


0.17 




^46 


15 


Al 


Mark 335 


0.45 


0.67 




^ 32 


16 


HI 


F 9 


1.05 


0.14 




^46 


30 


B1+ 


3C 249.1 


0.59 


0.17 




^ 32 


53 


B1++ 


3C 110^ 














'^ Equivalent width in A. 

" Lower limits to Lya/H/3 are estimated by the maximum contribution expected by a component of the same shift 

and width if peaking at 3(t the noise level. See text for a detailed explanation. 

'^ Consistent with intensity in all lines. 



Table 5 


. Very Broad 


Component (V 


BC): FWHM ~ 


10000 km 


s-i, A 


Vr ~ +2000 km s 


Sp. T. 


Name 






Intensity Ratio 






W" 




Civ/ 


Hen/ 


Cm/ 


Mgii/ 


Ha/ 


Lya/ 


Mgll/ 


Lya 






Lya 


Civ 


Civ/ 


Lya 


H/3 


H/3 


H/3 




A3 


I Zw 1 



















A2 


Mark 478 



















Al 


Mark 335 



















Bl 


F9 


0.36 


0.13 


0.18 


0.08 


1.26 


9.0 


0.68 


88 


B1+ 


3C 249.1 


0.67 


0.17 


0.17 






4.5 




41 


B1++ 


3C 110 


0.37 


0.10 




0.04 




5.0 


0.2 


27 




3C110^ 


0.50 


0.10 




0.07 




12.6 


0.83 


97 



Equivalent width in A. 

VBC only; see text for details. 
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Table 6. Interpercentilc Velocity Measurements 



Name 



Lya^" 



CivA1549=' 



A1900'' 



MgiiA2800'' 



H/3'' 



B 







R 



B 







R 







R 







R 



B 







R 



IZwl^ 4.2±0.7 15.5±3.4 1.1±0.9 1.7±0.8 2.4±1.3 ^0.4 

Mark478''3.4±l.l 33.0±2.4 1.3±1.2 1.6±0.6 5.1±1.0 ^0.4 

Mark 335*' 7.5±2.2 75.9±4.9 3.9±1.8 3.8±1.2 30.8±2.5 < 1.28 

Fairall 9^^ 4.0±1.5 28.5±3.2 7.2±2.6 2.2±0.7 16.3±2.4 2.7±1.0 

3C 249.1'^ 2.4±0.5 6. Oil. 2 1.5±0.7 1.4±0.2 3.6±0.6 1.1±0.3 

3C 110': 1.4±0.3 4.2±0.8 3.0±0.2 0.7±0.2 2.4±0.4 1.5±0.2 



0.2 0.8±0.3 


<0.4 


< 0.3 4.2±0.8 


0.4±0.3' 


0.3 0.6±0.4 


^0.4 


< 0.8 2.8±0.5 


<0.4 


0.9 6.1±1.8 


< 0.8 


< 0.7 4.5±0.7 


<0.5 


0.8 2.5±0.9 


< 0.5 


< 0.4 4.6±0.6 


0.5±0.4 



,0.2 0.8±0.2 0.2±0.1 



,0.1 0.5±0.1 O.lSitO.Ol 



0.15±0.10 2.4±0.3 ^0.1 
£0.1 3.2±0.2 <0.2 

0.15±0.10 4.0±0.2 0.2 ±0.1 
<0.1 2.2±0.2 0.6±0.3 

0.15±0.05 0.6±0.1 0.3±0.1 
O.litO.l 0.4 ±0.1 0.3±0.2 






"Fluxes in units of 10 ^^ergs s ^ cm '^. 
-5000 km s^l < Dr ^ -3000 km s^^; 0: 1000 km s^^ ^ Dr ^ 1000 km s^^; R: 3000 km s"! ^ i)r < 5000 km s'^. 
-7000 km s^l ^ i-r sg -5000 km s^^; 0: 1000 km s^^ ^ t)r ^ 1000 km s"^; R: 5000 km s^^ ^ t)r ^ 7000 km s"^ 
"^AllllAlSeO is measured for I Zwl and Mark 478; Cni]A1909 for the remaining objects. 

•: Strong Fen residual. 
®R: 3500 km s^^ ^ tir Sj 5500 km s~^ due to a blemish in the spectrum. 



Table 7. Black Hole Mass and Eddington Ratio 



Name 


fx 
at 5100 A 


log Lboi 


FWHM 


log 


Mbh 


logibo 

H/3 


iZ-'^Edd 




H/3 


H/3bc 


H/3 


H/3bc 


H/3bc 




[ergs s~^cm~^ J^^'^] 


[ergs s-i] 


[km s-l] 


[km s-i] 


[Mq] 


[M0] 






I Zw 1 


8.40E-15 


45.5 


1100 


1100 


7.3 


7.3 


-0.3 


-0.3 


Mark 478 


7.13E-15 


45.7 


1300 


1300 


7.5 


7.5 


-0.4 


-0.4 


Mark 335 


9.00E-15 


44.8 


1500 


15000 


7.2 


7.2 


-0.9 


-0.9 


Fairall 9 


7.52E-15 


45.3 


6500 


4550 


8.7 


8.4 


-2.0 


-1.7 


3C 249.1 


3.00E-15 


46.5 


8300 


5050 


9.5 


9.0 


-1.6 


-1.2 


3C 110 


1.50E-15 


46.8 


14000 


10700 


10.1 


9.9 


-1.9 


-1.6 
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Figure 1. Optical plane of the 4DE1 parameter space, i.e., FWHM{H/3) vs. -Rpdi. iJpcii measures the pr ominence of FellA45 70 with 
respect to H/3 (see text for definition). Occupied spectral bins are labeled. Grey dots are from the sample of lZamfir et al.l 1120081 ). Large 
black spots mark the position of the six sources considered in this study. 
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Figure 2. Multi-component fits of the Lya, ClvA1549, 1900 A, MgllA2800, H/3 spectral regions for the six object of this work ordered 
along the El sequence. All spectra are reduced to rest-frame and continuum subtracted. Ordinate is specific flux in units of 10~^^ ergs 
s"^ cm~^ A~^. The lower part of each panel shows residuals as a function of radial velocity. Thin black line: continuum subtracted 
rest-frame spectrum; thick dashed (magenta) line: multicomponent fit result. Thick continuous lines: blueshifted (blue), broad (black) 
and very broad (red) components. Thin orange line: narrow line components; thin dark green line: Fell template emission. In the 1900 
blend panel, the thick dark green line shows the adopted Felll template. In case of severe Felll contamination the Cin]A1909 line is not 
shown. 
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Figure 3. ClvA1549 profile (grey line), after continuum subtraction for the six sources considered in this study. The scaled profile of 
H/3 (H/3bc + H/3vBc) is shown as a thick line, and the residuals as a thin continuous line. For Fairall 9 and 3C249.1 it is possible to 
identify two major contributions to the residuals, namely one due to CivA1549nc, and one broad and shifted to the blue with respect to 
the source rest frame. The scaling is meant to show a maximum possible H/3 contribution to the ClvA1549 profiles; the actual ClvA1549 
retrieved from the multicomponent fits (see text) can be somewhat different from the one shown in this Figure. 



1550 



Rest Frame X[A] 
1600 1650 



1700 



1 1 1 r 




J I I I I I I I I I I I I I I I I I I I I I I t 



-10000 

AVj, [km s"'] 



10000 



Figure 4. The interpretation of the HellA1640 profile of 3C249.1. The "plateau" appearance of the far red wing of ClvA1549 is accounted 
for by a blueshifted (blue) and a redshifted very broad component (red) of HellA1640, whose shifts and widths match the ones of ClvA1549. 
The resulting profile is shown as the thickest black line. Units are as for Fig. [2] The dashed magenta lines traces the fit with all components 
added up. Other line styles have the same meaning of the previous Figure. 
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Figure 5. Predicted equivalent width of Lya in A (thick lines) and intensity ratio Lya/H/3 (thin lines) as a function of the ionization 
parameter U. The dashed lines refer to colunm density A'c= 10'^-'^ cni^^; the filled ones to Nc= 10^^ cm^^. 



APPENDIX A: NOTES ON INDIVIDUAL 
SOURCES 

I Zw 1 is the prototypical NLSyl source. Its UV spec- 
trum has been ana lyzed in detail by many workers, notably 
iLaor et alj (|l997bh . The MgllA2800 fit shows residuals that 
may be indicating an excess Fell emission around 2820 A. 

Mark 478 is spectroscopically rather similar to I Zw 1, 
with less extreme properties. The very low ClvA1549 BC 
intensity value may in part due to a narrow absorption close 
to the systemic velocity. 

Mark 335 belongs to bin Al. Its spectrum is consistent 
with the median spectrum of the bin, with low -Rpoii, weak 
blueshifted component. The profile of H/? is well described 
by a Lorentzian function, making this object similar to the 
rest of Pop. A even if -RfcII is as low as in Pop. B. 

Fairall 9 is a rather typical Bl source, with a relatively 
modest but appreciable VBC. The AhiiA1860/ Siiii]A1892 
and Siiii]A1892/CivA1549 ratios, along with _Rfcii«1, sug- 
gest some metal enrichment above solar. 

3C 249.1 is a lobe-dominated radio source. It shares many 
of the properties of several radio-loud quasars at higher red- 
shift, notably the strong CivA1549nc, and the shape of the 
Lya and of the ClvA1549-f-HellA1640 blend. 

3C 110 Some quasars show lines wi th strong VBC emis- 
sion that can be mistaken for a BC (ISulentic et al.ll2000d . 
e.g. PG1416-129). The correct interpretation for 3C110 may 
involve a single redshifted VBC component. The BC may 
be completely absent, as there is no Fell detection with 
the present data. Th is suggestion is also motivated by an 
analysis of 3C390.3 IjNegrete et al.l |2010|) which is a lobe- 
dominated RL quasar like 3C110. Physical conditions in- 
ferred for the BC of 3C 390.3 are similar to those of the 
VBC implying that there may be no BC in 3C 390.3. If 
the BC is completely suppressed in 3C110 then the entire 
ClvA1549 profile in that source may be VBC. A fit with a 
single, shifted Gaussian is worse than the one shown in Fig. 
(2] but not dramatically so. Then shift and FWHM values for 
a pure VBC ClvA1549 in 3C110 would be close to the canon- 
ical values (10^ and 10^ km s~^ respectively) for the VBC 
component in H/3. The spectrum of this object does lack a 
blueshifted component which might be closely associated to 
the BC. 



